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PREDICTION OF VIKING LANDER CAMERA IMAGE QUALITY 

Friedrich O. Huck, Ernest E. Burcher, Daniel J. Jobson, 
and Stephen D. Wall 
Langley Research Center 

SUMMARY 

Formulations are presented that permit predictions of image quality as a function 
of camera performance, surface radiance properties, and lighting and viewing geometry. 
Image quality predictions are intended to aid in diagnosir^ camera performance, in estab- 
lishing a preflight imaging strategy for the Viking lander mission to Mars, and in revis- 
ing this strategy if image data returned from Mars reveal unfavorable or unanticipated 
conditions. 

Predictions made for a wide range of surface radiance properties reveal that image 
quality depends strongly on proper camera dynamic range command and on favorable 
lighting and viewing geometry. Proper camera commands depend mostly on the surface 
albedo that will be encountered; to characterize the scene radiance initially, it is better 
to err with a dynamic range that is too wide to provide a sufficient number of quantized 
radiance levels for adequate pictures than with a dynamic range that is too narrow to 
encompass most of the scene radiance. Favorable lighting and viewing geometries depend 
mostly on the lander orientation with respect to the diurnal Sun path over the landing site, 
and tend to be independent of surface albedo and illumination scattering function. Side 
lighting with Sun elevation angles ranging from 10° to 30° is most often favorable for 
imaging spatial details and slopes. Within this range it is better to err with Sun elevation 
angles that may prove too high rather than too low for optimum imaging geometry. High 
Sun elevation angles are always favorable for measuring spectral reflectances. 

INTRODUCTION 

Two Viking spacecraft, each consisting of an orbiter and a soft lander, are scheduled 
to arrive at Mars in the summer of 1976 (refs. 1 and 2). Each lander carries two identical 
facsimile cameras together with many scientific instruments. 

The camera (ref. 3) features an array of 12 silicon photodiodes consisting of 4 broad- 
band channels with selectable focus for high-resolution imaging, 1 broadband channel for 
rapid surveys, 6 narrowband channels for multispectral imaging (color and near -infrared), 
and 1 narrowband channel for scanning the Sun. The instantaneous fields of view are 0.04° 
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for the four high-resolution channels and 0.12° for the other channels. The field of view 
ranges in elevation from 40° above to 60° below the horizon, and in azimuth to 342.5°. 

High sensitivity is obtained over a wide dynamic rai^e with only 6-bit encoding by the use 
of 6 commandable linear gains and 32 offsets. The camera scanning rates are synchro- 
nized with the lander data transmission rates of 16 000 bits per second to the two orbiters 
as relay stations and 250 bits per second directly to Earth. Image data can also be stored 
on a lander tape recorder. About lO"^ bits of image data will be transmitted during most 
days of the 60-day-long mission planned for each lander. 

It is generally known from Lunar Orbiter and Surveyor pictures of the Moon and 
from Mariner pictures of Mars and Mercury that good image quality depends not only on 
camera performance but also on lighting and viewing geometry and on surface radiance 
properties (i.e., albedo and illumination scattering function). Primary concerns for the 
Viking lander imaging investigation are degradation of camera responsivity by neutron 
radiation from the lander radioisotope-thermoelectric generators during the year-long 
cruise to Mars; uncertain lighting and viewing geometry despite efforts to control the 
lander orientation during final descent; uncertain surface relief and radiance properties; 
contamination of the camera window by dust adhesion and abrasion from dust storms with 
velocities up to 70 m/s; delay of several days between receipt of pictures and update of 
camera commands because of mission constraints; and initial trials and possible errors 
in selecting suitable camera gains and offsets. 

To aid the imaging investigation, this paper presents formulations that permit image 
quality predictions to be made as a function of camera performance, target radiance prop- 
erties, and lighting and viewing geometry. Image quality predictions as a function of ref- 
erence test chart features are intended to aid in diagnosing camera performance. Predic- 
tions as a function of surface features are intended to aid in establishing a preflight 
imaging strategy and in rapidly revising this strategy if image data returned from Mars 
reveal unfavorable or unanticipated conditions. This paper also presents results of pre- 
dictions as a function of lighting geometry for a wide range of surface radiance properties 
to aid in establishing broad guidelines for imaging. 

SYMBOLS AND NOTATION 


A area, m^ 

aQ,aj,a 2 physical parameters of illumination scattering function (see table V and 
eq. (48)) 

B(u,v,w) normalized image contrast of cone section (see eq. (42)) 


2 


b 


number of pixels contained in cone section 


c 

D 

e 

f 

G 

g 

h 

id 

ie 

JnO 

K 

k 

i^c 

i^CO 

kg 

^i 

ko 

L 

S. 


number of pixels per cone diameter 

lens diameter, m 

charge of an electron, C 

focal length, m 

channel gain 

phase angle, deg (see fig. 13) 
height, m 


photodiode noise-current density, A/yHz 


photosensor-array noise-current density, a/^1z 


nth order Bessel function 


photodiode signal current, A 
Boltzmann constant, J/K 

calibration factor for photosensor -array channels 
commandable offset interval, V (see fig. 7) 
gain constant, N^j^V 

ratio of signal integration time to sampling time 
fixed offset, V 

object distance from lens, m 
distance from lens in image space, m 
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N 


total rms noise, V 


digital number (see eq. (17)) 

Ng gain number (see fig. 7) 

Noff offset number (see fig. 7) 

Nq rms quantization noise, V 

N(X;t,e,g) object radiance, W/cm2-sr 
NER noise-equivalent radiance, W/cm2-sr 

ii surface normal 

o(x,i/) normalized spatial object radiance distribution 
R(X) spectral responsivity of photosensor array, A/W 

Rf preamplifier feedback resistance, 

Rj channel resistance, 

r radius, m 

S(X) solar irradiance, W/cm2-/j,m 

S(x,i/^) camera signal, A 

s dimensionless variable for spatial frequency (see eq. (4)) 

T temperature, K 

u dimensionless variable for defocus (see eq. (6)) 

V photosensor -array output voltage, V 

V dimensionless aperture radius (see eq. (5)) 
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W noise-equivalent bandwidth, Hz 

w dimensionless variable for signal integration time (see eq. (lib)) 

X,Y camera elevation and azimuth sampling interval, respectively, rad or deg 

a surface or cone slope, deg 

/3 instantaneous field of view, rad or deg 

y(X) mirror reflectance 

6 increment of cone azimuth angle, deg 

e emittance angle, deg (see fig. 13) 

? azimuth angle between incident light and cone surface element normal, deg 

(see fig. 13) 

azimuth angle defining cone shadow limits, deg (see fig. 14) 

6 azimuth angle between incident light and camera viewing axis, deg 

L incidence angle, deg 

K number of quantization levels 

A tribar width on reference test chart, m 

X wavelength, fim 

v,(jj spatial frequency along elevation and azimuth direction, respectively, Ip/rad 

p(X) albedo 

t(x) spectral transmittance 

point- spread function 
t( ) spatial frequency response 
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optical thickness of atmosphere 
illumination scattering function 
mirror scan rate, deg/s 

azimuth angle between cone surface normal and camera viewing axis, deg 
(see fig. 14) 

azimuth angle defining limits of visual shadow, deg (see fig. 14) 
angle of numerical aperture, deg (see fig. 8) 
sampling or comb function 
rectangular function 

Subscripts: 

a atmosphere 

c camera 

cn cone 

cw camera contamination window 

e electronic 

f spectral filter 

i image 

S. lens 

m Mars 

0 flat surface 

p photosensor aperture 


T 

4>U,e,g) 

X 


^sh 

mo 

no 
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pw photosensor -array window 

sh shadow 

V visible 

w camera window 

+ brighter than flat surface 

darker than flat surface 
Abbreviations: 

BB broadband 

IR infrared 

Ip line pair 

SEA Sun elevation angle 


The symbol A in front of a parameter indicates an increment of that parameter. 

The bracket ( ) around a parameter indicates that parameter which is to be 
estimated. 

A bar over a symbol denotes an average or weighted value. 

A circumflex ('') over a symbol denotes a Fourier transform. 

CAMERA 

Description 

Location .- Figure 1 presents a view of the Viking lander. Its two cameras are sep- 
arated by 0.8 m and view the surrounding terrain from a height of 1.3 m. Figure 2 shows 
an outline drawing of the field of view available to the cameras. A large segment of the 
camera field of view is obscured by the lander structure. Important to calibrating the 
camera performance are three reference test charts, of which each camera can view two. 
The shaded surface area defines the area accessible to a surface sampler for biological 
and organic and inorganic chemical analysis. 
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A roll maneuver (i.e., rotation about an axis parallel to the camera axis) is planned 
for the final lander descent phase. An important constraint on lander orientation is to 
avoid shadowing the sampler area at favorable Sun elevations. A general goal is to obtain 
morning and afternoon lighting of the sampler area, first from one direction and then from 
the other, to insure that all features in this area can be viewed at one time or another out 
of shadow. However, contrary to the best efforts of the flight engineers, the spacecraft 
may land in an orientation that results in a different, less favorable lighting geometry. 

Construction . - Figure 3 presents a simplified cutaway view of the camera. The 
upper elevation assembly rotates in azimuth during operation, and the lower azimuth 
assembly is bolted to the lander. General design and operating characteristics are listed 
in table I and performance characteristics in table II. 

The thin camera housing and the major support structure are made of beryllium to 
provide lightness, stiffness, and low thermal expansion and contraction. The inner sur- 
faces of the housing are treated with a matte black finish to reduce internal light reflec- 
tions, and the outer mating surfaces are gold plated to suppress electromagnetic interfer- 
ence. A 3- to 5-cm-thick blanket of fiberglass and aluminum foil insulates the entire 
camera to reduce radiative and convective transfer to the cold Martian atmosphere. In 
addition, thermostatically controlled heaters are located in the vicinity of critical compo- 
nents to maintain their temperature above -40° C. 

The elevation assembly is hermetically sealed and filled with argon gas at a pres- 
sure of nearly 1 bar (105 N/m2). The azimuth assembly is exposed to the Martian atmo- 
sphere. However, dust is prevented from entering by a labyrinth seal between the rotating 
and stationary parts of the camera. 

The camera has two 2.5-mm-thick fused silica windows, each with an outer coating 
of Si02, which discourages dust adhesion and abrasion, and an inner coating of MgF2, 
which reduces internal light reflections. The outer window, referred to as the contami- 
nation cover and not shown in figure 3, is hinged on one side and spring-loaded. By rotat- 
ing the camera nearly 360° in a rapid slew mode, a lever can be forced against the side of 
the dust post with sufficient impact to release a catch on the contamination cover; a spring 
then moves the cover out of the optical path. This cover protects the inner window from 
gaseous condensation durir^ the lander heat sterilization process prior to flight, and from ^ 

dust contamination during the first few days or weeks of operation on Mars. | 

I 

The narrow window entrance usually hides behind a post to avoid unnecessary expo- 
sure to dust. The post contains two small tungsten filament light sources that permit the 
camera performance to be checked when the lander is stored inside the Viking spacecraft 
during the flight to Mars. Carbon dioxide carried under pressure in a container inside the i 
lander can be released on command through a tube, which is attached to the post, to blow 
dust off the window. 

! 

II 
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Light which passes through the window is reflected by a mirror that nods up and 
down to provide the elevation scanning; scan profile and rates are given in figure 4. 

Image data are acquired only during the relatively slow upward motion. The reflected 
light is imaged by an achromatic triplet lens onto a photosensor array. The lens aperture 
diameter is 0.95 cm and the focal length is 5.37 cm. light falling on the photosensor 
array is transduced into an electrical signal. 

Immediately below the lens is a dark reference shutter which consists of a rotary 
solenoid that controls a paddle-shaped black metal shutter. When the solenoid is ener- 
gized, the solenoid shaft moves the shutter over the exit pupil of the lens to block external 
light from the photosensor array. The design includes a spring that forces the shutter to 
the open position in case of a failure of the solenoid or associated electronics. The dark 
reference shutter is closed for sampling the photosensor-preamplifier dark current, which 
is stored for later subtraction from the signal current, and for calibrating the photosensor 
array with a small tungsten filament lamp located below the shutter. The inner surface 
of the optical tunnel between the lens and the photosensor array contains a multiaperture 
baffle that is designed not only to reduce veiling glare but also to attenuate radio -frequency 
interference from the lander antennas. 

The elevation assembly, which contains the camera elements described so far, tap- 
ers down to a narrow shaft below the photosensor array. All the cables that provide con- 
nections between the electrical elements in the elevation assembly and the control and 
signal processing electronics in the azimuth assembly spiral down through the inside of 
this shaft. The azimuth servo assembly, which is mounted on the outside of the shaft, is 
similar to the elevation servo assembly mounted on the mirror shaft. Both assemblies 
consist of a torque motor, a tachometer, and a resolver. 

Signal flow .- The light-sensitive part of the photosensor array is a miniature array 
of 12 silicon photodiodes. Their responsivities are plotted in figure 5. Four of the diodes 
are spaced at different distances from the lens for electrical focus selection for high- 
resolution (black-and-white) imaging. Six of the other diodes are covered with optical 
filters for multispectral (color or infrared) imaging. Another diode is unfiltered and 
used for rapid surveys of the terrain, and one is covered with a red filter and used for 
viewing the Sun directly. The nominal aperture diameter is 0.041 mm for the high- 
resolution diodes, and 0.119 mm for the low-resolution diodes. 

The cameras are exposed to a steady flow of neutrons and gamma rays from the 
lander radioisotope-thermoelectric generator during the nearly year-long cruise from 
Earth to Mars. Radiation damage by the neutrons degrades the photodiode performance, 
particularly in the infrared region from 0.8 to 1.1 jum, where the absolute responsivity 
will be reduced by an estimated 15 to 40 percent and the temperature dependence of the 
relative responsivity will be increased. 
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Figure 6 shows a simplified circuit diagram of the photosensor array and video 
processing electronics. Eleven of the diodes have their own preamplifier to convert their 
small output current into a usable signal, while the Sun diode signal current is sufficiently 
high to be directly handled by the video processing electronics. The preamplifiers are 
conventional current-to -voltage converters. Their outputs are summed into a buffer 
amplifier. The gain of this stage is the ratio of the amplifier feedback resistance to the 
input resistance of each channel. The input resistances were selected to compensate for 
different photodiode aperture sizes and filter transmittance characteristics. 

The output from the summing amplifier is passed to circuits which provide com- 
mandable gains and offsets, automatic dark-current subtraction, and analog-to-digital 
(A/D) conversion. The primary electrical filtering which the signal undergoes is the 
integration in the A/D converter. 

Six gains can be commanded in steps of powers of 2, and 32 offsets can be com- 
manded in small steps from a slightly negative offset up to nearly the full dynamic range 
of the camera, as shown in figure 7. These commandable gains and offsets provide a 
means of decreasing quantization intervals (at the expense of dynamic range) without 
increasing data bits per pixel. 

Control logic.- Camera operation is controlled by a 48-bit command word from the 
lander guidance, control, and sequencing computer. When the camera logic senses that 
the lander computer has turned on power for the camera, it clears all its memory ele- 
ments and commands the azimuth and elevation servos to a preset position. Immediately 
after, the azimuth start and elevation pointing angles are loaded into their respective 
counters, forcing the servos to slew to the commanded positions. Upon receipt of clock 
pulses from the lander data acquisition and processing unit, the camera logic selects first 
engineering data (such as temperature) from various sources and then video data. Video 
data from the A/D converter are passed through a 524-bit memory which converts the 
75-percent mirror scanning efficiency to a 90-percent video data format efficiency. This 
process is repeated for each line scan in an image. Three imaging modes can be 
commanded: 

(1) A high- resolution imaging mode is effected by sending one incrementing pulse to 
the elevation servo between video samples and to the azimuth servo between active line 
scans. 

(2) A low-resolution (survey) imaging mode is effected by tripling the number of 
incrementing pulses sent to the elevation servo between video samples and to the azimuth 
servo between active line scans. 

(3) A multispectral (color or infrared) imaging mode is effected by alternately 
selecting three diodes (either blue -green-red-blue, etc. or IR3-IR2-IR1-IR3, etc.) for 
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successive elevation scans and inhibiting the azimuth servo until three elevation scans 
have been completed. 

In addition, repeated line scans are effected by inhibiting the azimuth servo. This 
mode is used for radiometric measurements of atmospheric and surface features, and for 
observing variable features. 

Data are acquired at either one of two rates: synchronous with the lander data 
transmission rate of 16 000 bits per second to one of the two Viking orbiters, or 250 bits 
per second directly to Earth. Data acquired at either rate can also be stored on a lander 
tape recorder, permitting the cameras to be used during favorable imaging opportunities 
independent of data transmission periods. Typical picture formats and their data content 
are listed in table HI. 


Analytical Model 

Imaging process .- The process by which the facsimile camera transfers the (con- 
tinuous) object radiance distribution o(x,^) into a (discrete) electrical signal S(x,i/^) 
can be approximately formulated by the expression (ref. 4) 




( 1 ) 


where K is the photodiode current for uniform radiance; the camera point- 

spread function; x nnd xf/, the angular elevation and azimuth imaging coordinates, 
respectively; and X and Y, the angular elevation and azimuth sampling intervals. 


respectively. 


The symbol ♦ 


denotes convolution, and 




denotes the sampling 


(ref. 5) or comb (ref. 6) function. The azimuth sampling interval Y is equal to the cam- 
era azimuth stepping interval times cos x» where x is measured from a plane normal 
to the optical axis of the objective lens. 


An approximation is introduced into the formulation of equation (1) by the separation 
of spatial and spectral object and camera characteristics. Actually, o(x,>/^) and Tg(x,i^) 
are functions of wavelei^th, and the spatial convolution should therefore be integrated over 
wavelength. However, it is generally convenient to separate spatial and spectral charac- 
teristics, and let the average photosensor signal current K account for the spectral 
characteristics. This approximation permits o(x,i^) and Tj,(x,i//) to be expressed as 
normalized functions, whereas S{x,x}/) takes on the dimension of K. It requires that all 
spatial computations be performed at the surface-radiance and camera-responsivity 
weighted wavelength X. 
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The optical-mechanical line-scan imaging process of the facsimile camera is 
implicitly a function of time. The formulation of equation (1) implies therefore that the 
convolution of the object radiance distribution with the camera point-spread function be 
performed as a function of the time each pixel is acquired in a picture (ref. 7). How- 
ever, if neither object radiance distribution nor camera response varies with time (as is 
assumed here), then it is immaterial whether the pixels in a picture are formed simul- 
taneously or sequentially. 

The imaging process of the facsimile camera is more conveniently evaluated in the 
frequency rather than spatial domain. The Fourier transform of equation (1) can be 
written as 


S(i^,co) = K o(i',co) T^(v,u}) * XY rn(Xi^,Yw) 


( 2 ) 


where v and w are the spatial frequency along the elevation and azimuth direction, 
respectively, and the circumflexes represent the respective Fourier transform variables 
of the functions in equation (1). After this signal is passed through a proper low-pass 
filter, 


S(i^,a;) » K o(i^,U)) Tg(i^,w) 


(3) 


The approximation given by equation (3) disregards all frequency components of the side- 
bands that may fall into the passband of the filter and cause image degradation due to 
aliasir^. (See refs. 3, 8, and 9 for details.) 

Spatial response .- The spatial frequency response of the camera is the product 
of the frequency responses of the objective lens t^, photosensor aperture fp, and signal 
electronics Tg. To formulate and evaluate the camera frequency response, it is conve- 
nient to change from spatial variables {v and a>) to dimensionless variables (s, v, 
and u). The latter take advantage of the circular symmetry of the objective lens and 
photosensor aperture, and allow numerical results of the following equations to be readily 
used for a variety of specific configurations. 


Spatial frequency is accounted for by the dimensionless variable s 
(ref. 10) 


given by 


sin Q, 


D 


(4) 


where X is the surface-radiance and camera- responsivity weighted wavelength; the 

distance from the image plane to the lens; sin the numerical aperture; and D, the 
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diameter of the objective lens. (See fig. 8.) The photosensor aperture radius rp and 
the instantaneous field of view ^ formed by this aperture are accounted for by the 
dimensionless variable v given by (ref. 10) 


V = is rp sin 

X ^ x\2fpyP l\2)^ 


(5) 


where £p is the distance between the lens and the aperture. The amount of defocus 
introduced by differences in image-plane and aperture-plane distances from the lens is 
accounted for by the dimensionless variable u given by (ref. 10) 


u = 


is 

X 


Af tan sin ® — A£ 

2X 



( 6 ) 


where 



Object-space distances L are related to image-space distances 


£ by the thin-lens formula 


i=i+i 

f L £ 


(7) 


where f is the lens focal length. 

The lens performance is essentially diffraction limited over the angular extent of 
the photodiode array. The frequency response of a defocused, diffraction-limited lens 
has been formulated by Hopkins (refs. 11 and 12) as 


T^(u;s) = ^ us2^|yJl(us) + ^ sin 2y[ji(us) - J3(us)j - sin 4y|j3(us) 

- J5(us)j + . . .^ - ~ sin^i us2^|sin y|jo(us) - J2(us)j - i sin 3y|j2(us) 

- J4(us)j + i sin 5yjj4(us) - J0(us)j - . . ^ 


(8a) 


— 1 s 

where y = cos ^ . As the photosensor aperture approaches focus (u = 0), the lens 

response becomes 
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Tj^(0,s) = |(2y - sin 2y) 


(8b) 


Table n lists the in-focus object distances Lp for the photodiodes. Numerical results 
of equations (8) for various amounts of defocus are presented in figure 9(a). If the best- 
focused high-resolution diode (see table II) is used at each distance, then u = 15 for 
object distances of 1,7 m to infinity; for the low-resolution diodes, u = 45. 

The photosensor apertures are circular, and their frequency response is given by 


Tp(vs) = 


2Jj^(vs) 

vs 


(9) 


Numerical results of equation (9) are plotted in figure 9(b). 

The primary electrical filtering occurs in the running mean integrator of the A/D 
converter. The frequency response is 


Tp = sine 


(xTii^) = 


= sine kiX!^ 


(10) 


where sine x = the mirror line-scan rates (given in fig. 4); Ti, the integra- 

7TX 

tion time; and kj^ (=0.44), the fraction of the sampling interval X = ^ during which the 
signal is integrated. Substituting equation (4) for v (with cj = 0) yields 


Tg(ws;ct)=0) = sine ws 


(11a) 


where 

w = — kt/3 (11b) 

2X 

Numerical results of equations (11) are plotted in figure 9(c). 

The spatial frequency response of the camera as a function of the dimensionless 
variables is 

T^.(u,v,w;s) = T^(u,s) Tp(vs) Tg(ws;w=0) (12) 

Numerical results of equation (12) are plotted in figure 9(d). 
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Spectral response .- The photodiode current K(t,e,g) for a spatially uniform radi- 
ance N(X;t,e,g) is predicted by (ref. 10) 

K(t,e,g) N(X;t,e,g) t^CX) R(X) dX (13) 

where t^.(x) is the optical throughput of the camera, and R(X) and kg are the 
responsivlty and calibration factor, respectively, of each photosensor-array channel. 

The radiance of a surface is 

N(X;L,e,g) = I S(X) Tg^(io;^') pW 0(t-,e,g) (14) 

where S(X) is the solar irradiance; '^a(‘-o>^)’ atmospheric transmittance; and p(X) 
and <p{L,e,g), the surface albedo and illumination scattering function, respectively. 

The optical throughput of the camera is 


Tc(X) 


cw 


(X) T^(x) y(X) T^(x) 


y(X) T^(X) 


(cover in place) 
(cover removed) 


(15) 


where transmittance of the contamination cover, 

window, and lens, respectively, and y(X) is the reflectance of the mirror. Values are 
plotted in figure 10. 

The responsivity for each photosensor-array channel is 


R(X) - 'Tp-w(X) T£(X) R (X) 

where Tp^(x) is the transmittance of the photosensor-array window (about 0.96); t^(x), 
the spectral filter (if present) transmittance; and R'(X), the photodiode responsivity. 
Values of R(X) are plotted in figure 5. 

The photosensor output voltage V(t,e,g) for a photodiode current K(t,e,g) (see 
fig. 6) is 


V(t,e,g) = K(t,e,g) RfG 

where Rf is the preamplifier feedback resistance and G the channel gain. 


(16) 
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The calibration factor k^, for each camera channel is selected so that predicted 
photosensor-array output voltages V, from equations (13) and (16), agree with camera 
calibration measurements. Camera response and calibration data are given in 
appendix A. 

The quantized signal, or digital number Njj, is related to an analog signal V by 
(see fig. 7) 


Nd = ^(V - kcoNoH + ko) 


(17) 


where Ng is the commandable gain number (ranging from 0 to 5); Ngff, the command- 
able offset number (ranging from 0 to 31); lig, the gain constant (in 
offset interval; and ko, the fixed offset. 

The surface-radiance and camera- responsivity weighted wavelei^th X at which 
all spatial computations should be performed is given by 


X = 


pOO 

X S(X) T^(to;>) p(X) R(^) 

^oo 

J S(x) Ta(to;x) p(x) t^(X) R(X) dX 


(18) 


Noise .- The two most significant electronic noise sources in the photosensor array 
are shot noise in the photodiode current and Johnson noise in the preamplifier feedback 
resistor. Their equivalent noise-current density (in at the input of the preampli- 

fier is (ref. 3) 


i 


e “ 




(19) 


where e is the charge of an electron; i^j, the photodiode current; k, the Boltzmann 
constant; and T, the absolute temperature. It is generally complicated to account rig- 
orously for variations in shot noise as a function of variations in the photodiode current; 
instead, an average value for the shot noise based on the signal current K = K(t=e=g=0) 
(i.e., for 0(t,e,g) = 1) is used as a conservative approximation. Figure 11 presents a 
plot of the total photosensor noise -current density ig as a function of the signal cur- 
rent K. 

The quantization of the electrical signal for transmission is a basic limitation of 
digital data in providing the true value of a signal, just as random noise is a limitation of 
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analog data. For k(= 64) uniformly spaced quantization levels within the dynamic range 


AV = 




kg 


2 the effective mean-square quantization noise (ref. 13) is 





( 20 ) 


When electronic noise and quantization noise are combined, the total rms noise 
magnitude referred to the photosensor -array output becomes 


N=y(ieRfG)^W + Nq2 (21) 

where W is the noise -equivalent bandwidth of the video processing electronics 
(W = 2.8 kHz for the rapid scan rate and W = 55 Hz for the slow scan rate). 

Sensitivity .- A common measure of sensitivity is the signal-to-rms-noise ratio 
V/N, where V is the photosensor -array signal voltage given by equation (16) for 
(f>{L,e,g) = 1, and N is the rms noise given by equation (21). A related performance 
parameter is the noise -equivalent radiance NER given by 



where N(X) is given by equation (14). Typical values are given in table IV. 

Diagnosis 

If the quality of the initial pictures received from the lander is adequate and roughly 
concurs with preflight predictions, then the basic preflight imaging strategy can be con- 
tinued with only minor changes to optimize camera commands for subsequent investiga- 
tions and to respond to terrain features of special interest. However, if the quality of the 
pictures is poor, then the preflight imaging strategy may need to be revised. Poor image 
quality could arise because of unanticipated surface properties, unfavorable lander orien- 
tation, or degraded camera performance. 

A computer program for image quality prediction should aid first as a diagnostic 
tool in isolating the cause of picture degradation, and thereafter as a predictive tool in 
revising the imaging strategy. This diagnostic and predictive fimction can be accom- 
plished as follows: If image data of a reference test chart agree with predictions, then 
the camera must be operating properly and the poor quality of pictures of the terrain 
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must have been caused by unanticipated surface radiance properties. However, if image 
data of a reference test chart disagree with predictions, then the poor quality of all pic- 
tures will probably have been caused, at least in part, by a degradation in camera per- 
formance, Whatever the initial pictures and camera engineering data may reveal, the 
resultant conjectures about degraded camera performance and unanticipated scene radi- 
ance properties can be accounted for with a prediction computer program by altering 
pertinent camera and scene parameters until predictions come into agreement with the 
image data. Thereafter, the prediction program can again aid in establishing favorable 
camera commands and lighting and viewing geometries. 

To accomplish the diagnostic function, it is necessary to predict camera perfor- 
mance in terms of reference test chart features. Three identical reference test charts 
are located on top of the lander for this reason. Each camera can view two of these 
charts at a distance of about 1 m. A reference test chart is shown in figure 12; it pro- 
vides 11 reflectance references (gray patches) for radiometric calibrations, 3 color 
patches to aid reconstruction of color images, and 3 tribars to check the camera fre- 
quency response. 

The reflectances of the gray and color patches are nearly Lambertian for light 
incidence angles ranging from 20® to 60®; the reflectances of the gray patches are also 
nearly constant, within ±5 percent, with wavelength from 0.4 to 1.0 jitm. The absolute 
reflectances given in figure 12 can be used with these simplifications with a 2 -sigma 
accuracy of ±9 percent. Absolute radiometric camera performance predictions may be 
checked by substituting cost for 0(t,e,g) and absolute reflectance values for p(X) 
in the surface radiance equation (14), and corrected by changing camera parameters until 
predictions and camera data agree. Camera parameters which may typically be adjusted 
are the window transmittance t^(X) or the photosensor-array calibration constants k^. 

The reflectances of the gray patches relative to each other are accurate to within 
±4 percent for any incidence angle and wavelength. This higher accuracy is important 
because the ratio of reflectances will noticeably change with appreciable dust adhesion or 
abrasion and will thus give warning to be careful about the use of the reference test chart 
as an absolute reflectance reference. 


The three tribars of the reference test chart provide an aid in checking the camera 
frequency response. The camera signal is predicted as a function of its (one -dimensional) 
sine-wave frequency response ^^(i.') along the elevation line-scan direction by the 
expression 


AV = ^ V(t) 


o 

I 

i=l 

j=2i-l 


Jbll 


i+1 


Tj.(u,v,w;js) 


sine 


/3L 

57.3A 


(23) 
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Sufficient accuracy is obtained with the first three terms of the expansion. The sampling 
process generates a signal which is statistical rather than deterministic in nature. The 

factor sine accounts for an average reduction in signal contrast; the actual signal 

57. 3A 

may have a slightly higher or lower contrast. The tribar widths A are 12.8, 6.5, and 
4.4 mm per line pair (corresponding at 1 m to 1.4, 2.7, and 3.9 Ip/deg, respectively). 
Appendix B presents a comparison of predicted camera square-wave responses with 
measurements. 


IMAGE QUALITY CRITERIA 
Definition 

A single criterion cannot completely describe the quality of an image; the most that 
one can hope to do is to define a criterion that is satisfactory for a particular objective. 
The objective of this study is to optimize lighting and viewir^ geometry and camera gains. 
Within this constraint it can be said that it is desirable to resolve small spatial details 
and slope variations in the broadband (survey and high-resolution) imaging modes, and 
subtle spectral reflectance variations in the narrowband (color and IR) imaging modes. 

The capability to resolve small spatial details and slopes is defined here as the 
minimum detectable cone diameter and cone slope with respect to a level surface. An 
upright cone with surface properties of the surrounding terrain seems intuitively repre- 
sentative of many features and has no preferred surface orientation azimuthally about its 
axis. If the cone angle is chosen to be steep, a condition yielding high surface contrast, 
then the detectability of this target becomes primarily a measure of the camera capability 
to resolve small detail. If the cone ar^le is chosen to be shallow, a condition yielding low 
surface contrast, then the detectability of the same target becomes primarily a measure 
of the camera capability to resolve shallow slopes. Cones have previously been used for 
image quality predictions for the Ranger (ref. 14) and Lunar Orbiter (ref. 15) missions. 

The capability to resolve surface reflectance variations in each one of the camera 
spectral channels is defined as minimum detectable albedo difference, that is, that differ- 
ence in albedo which results in a given signal-to-noise ratio, a level surface being 
assumed. 


Formulation 

Lighting and viewing geometry .- Figure 13 defines the angular relationships between 
surface slopes and lighting and viewii^ geometry. The illumination scattering function 
<^(t,€,g) accounts for the dependence of surface reflectance on this geometry as a func- 
tion of the ai^le t between incident radiation and surface normal, the angle e between 
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emitted radiation and surface normal, and the phase angle g between incident and emit- 
ted radiation. These three angles can be determined from the following relationships: 

cos g = cos 6o cos Lq + sin €q sin Lq cos 6 (24) 

cos L = cos a cos Lq - sin a sin Lq cos ? (25) 

cos e = cos a cos eQ - sin a sin 6o cos (26) 

where Cq and Lq are the emittance and incidence angles, respectively, for a level sur- 
face ^with normal ho), a is the angle of a surface element with respect to a level sur- 
face (and with normal h), 0 is the azimuth angle between incident radiation and the 

camera optical axis, and ? and are azimuth angles of the plane formed by the nor- 
mals h and ho with respect to incident radiation and the camera, respectively. The 
three azimuth angles are related by 

0 + ? + ^ = 271 (27) 

Cone geometry .- This formulation generally follows the geometric relationships 
derived by Rindfleisch and Willingham (ref. 14). However, their relationships were con- 
strained to nearly vertical viewing geometries so that no part of the cone, including its 
shadow (if present), could be obscured from the camera view by the cone itself. This 
constraint is not valid for a camera viewing a scene from a lander. Relationships from 
reference 14 are therefore generalized to account for those viewing geometries in which 
part of the cone or its shadow is obscured from the camera view by the cone itself. 

Figure 14 illustrates the geometry of an upright cone with a base of radius r^j^, 
a height hcn> a slope a defined so that 

a = tan"l— fo < oi < I'j (28) 

’^cn \ 

h is given in terms of the azimuth angle increment 6 

(29) 

If Lq = 90° - Q!, a shadow is cast by the cone on a level surface (as shown shaded 

in fig. 14(a)) for -Ci, < ? < Ccv,, where C. = cos"^ ^ ; this shadow has an area 

tan a tan Lq 


A surface element with normal 
as 

.2 


AA = 


5r 


cn 


2 cos a 
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(30a) 


•^sh “ ^sh ■ ?sh) 

Similarly, if = 90° - a, a "visual" shadow is formed (fig. 14(b)) which obscures 
part of the cone and possibly part of the cone shadow (if there is one) from the camera. 

It is assumed that the rays of light intercepted by the camera are parallel instead of con- 
vergent, since only small targets are of interest and the convergence angle is therefore 
small. Under this condition, area elements AA of the cone not visible to the camera 
are defined by the range ~4^sh < ^ < ’^sh’ ''^here = cos"^ . The shadow 

u3.ll Of U3I1 6q 

area still visible to the camera becomes (see fig. 14(c)) 


■^v,sh 



(30b) 


where 4 = |(Csh + ’^sh ' ®)‘ 

Several geometries must be accounted for when equation (30b) is used. The azimuth 
angle 6 between incident radiation and the camera optical axis must always be chosen 
so that 6 < TT. When 6 = + i/^gjj, the "visual" shadow does not overlap any part of the 

illumination shadow, and equation (30b) reduces to equation (30a). When ^ ^sh ’^sh’ 

the illumination shadow is completely contained within the visual shadow, and the angle ^ 
must be set equal to ?gjj so that Agj^ = 0. When 6 + = ?gj^, the visual shadow is 

completely contained within the illumination shadow, and the angle ^ must be set equal 
to l//gh- 

Target and background are considered to be of the same material, so that the vari- 
ation in cone reflectance with azimuth angle is dependent only on the illumination scatter- 
ing function. But it would still be exceedingly difficult to translate this reflectance vari- 
ation exactly into an image signal. The reflectance variation is therefore approximated 
by an average reflectance of that part of the cone which has a higher -than-background 
reflectance and another part which has a lower -than-background (including shadow) 
reflectance. 

The total projected area of all cone sections which have a higher-than-bacl^round 
reflectance and are visible to the camera may be expressed as 


= AAYcos6i=-^ii^ycose 

4^ 2 cos a. Lu 




(31) 


where The average value of the illumination scattering function of this 

area becomes 
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(0i(t ,€,g) > 0o) 


(32) 


<^+ = ^^0i(t,e,g) cos ej 
i 


J'i>i(‘.e,g) cos «i 


y 


cos 




Similarly, the total projected area of all cone sections, including the shadow cast by 
the cone on a level surface (if present), which have a lower -than-bacl^round reflectance 
and are visible to the camera may be expressed as 


A_ = AA ^ cos ej + Agj^ cos = r 


2 

cn 



tan 


^sh 


^sh - tan 4 + cos 


(0^(t,€,g) < 0o] (33) 


The average value of the illumination scattering function of this area becomes 


6 ^0.(t,€,g) cos 

cose^=-^ 

f 

(<^i(t,e,g) < <Po) (34) 

Reflectance inside the shadow is assumed to be negligible. 

Cone signal. - Several definitions and assumptions are made to simplify 
formulations: 

(1) According to equation (14), the maximum surface radiance for a given albedo 
(i.e., for 0(t,e,g) = 1) can be defined as 


«o(^ 


COS Sil +2 cos a cos €^(tan - tan | + ^ 


N(X) = i S(X) Ta(X) p(X) (35) 

The corresponding maximum photosensor signal current and photosensor-array output 
voltage are defined, respectively, as K and V. 

(2) The radiance variation of the cone surface can then be defined as 

ANj^(X) = N(X) (36) 
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and the corresponding variation in photosensor-array output voltage as 


AV^ = V = KRG 


( 37 ) 


(3) The two cone areas A . and A_ are assumed to be circular in shape with 
radius r^ and r_, respectively; that is, A^ = Trr^ . Computational requirements are 
thus reduced, since circular symmetry permits the use of the one -dimensional Hankel 
transform instead of the general two-dimensional Fourier transform. This assumption 
is justifiable because most of the radiance contained in the spatial frequency distribution 
of the cone radiance is generally clustered around the lower spatial frequencies. The 
assumption is also consistent with the earlier approximation of cone radiance variations 
by two average values. 

(4) The cone radius is expressed in terms of the camera instantaneous field of view 
/3 so that r^jj^ = - /3L, where L is the distance from target to camera. The product 

|8L is the diameter of a pixel, and c is the number of pixels per cone diameter. The 
diameters b^. of the areas A^. can thus be expressed as a function of pixels as 






(38) 


where the areas A^. are given by equations (31) and (33). 

From these definitions and assumptions, the spatial frequency distribution of the 
cone radiance (in image space) is given by the Hankel transform of the two areas A^., 
assumed to be circular with radius b^^v, as 


AN(x;b±vs) = AN±{X) J --- 
' pzn pb 

Jq *^0 


pb.v 

2ti \ Jq(v's) v' dv' 


pb.v 

' v' dv’ dw 


^ b^vs 


(39) 


where s and v are given by equations (4) and (5), respectively. 

The effect of the camera frequency response in reducing the contrast of spatial 
detail is accounted for by the inverse Hankel transform of the product of the cone radi- 
ance frequency distribution AN^X;bj_vs) and the camera frequency response Tg(u,v,w;s), 
as given by 
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/ N 2J,(b.vs) ^ 

AV±(u,v,w;b±j = AV± ' ^ys 


s ds 


(40) 


Upon averaging, as a final approximation, over the circular image area of radius bj.v, 
equation (40) becomes 


pb±V p«> ZJjID.VSJ 

AV^^u,v,w;b±j = AVj. J j — ^ — - Jo(v's) s ds v' dv’ 


2Ji(bj.vs) 


bj.vs 


Equation (41a) can be rewritten (ref. 16) as 


(41a) 



B , (u,v,w) = 


b, ^v2 


pOO 


2J 




b^vs 


T„(u,v,w;s) s ds 


(41b) 


(42) 


Mnimum detec table cone diameter and slope.- The detectability of a cone depends 
on its image contrast and size and on noise. The contrast is given by the foregoing 
results as 


AVj.(u,v,w) = AV^ B(u,v,w) = V A<f)^ B(u,v,w) (43) 

The size is given by two circular areas with radii b±v. According to the central-limit 
theorem, the mean of a sample of size n has a standard deviation Oy/|(n if o is the 
population deviation. Here a is the total rms camera noise N given by equation (21), 
and the number of independent samples n is the ratio of image area to pixel area; that 
is. 


I 
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( 44 ) 



The effects of contrast, size, and noise on the detectability of a cone can be com- 
bined in an expression for signal-to-noise ratio as 


{ ^ ^[^V+(u,v,w) b+ + AV_(u,v,w) b_j (45a) 

{n/^ ^ B+(u,v,w) + A0.b. B_(u,v,w)J (45b) 

If it is desirable to determine the minimum detectable cone slope, the procedure should 
be to specify c and then to find the slope O! which yields a given value of (S/N). Or 
if it is desirable to determine the minimum detectable cone diameter, the procedure 
should be to specify a and then find the normalized cone diameter c (in pixels) which 
yields a given value of (S/N). Both procedures require iterative evaluations of equa- 
tions (45) until a solution with desired accuracy is found. The signal-to-noise ratio (S/N) 
is somewhat arbitrarily taken to be 10. 

Minimum detectabl e albedo difference .- Three assumptions are made: (1) The sur- 
face is assumed to be level so that <^(t,e,g) = surface albedo is assumed 

to be constant over the narrow wavelength range of the camera spectral filters; and (3) an 
albedo difference can be detected for a signal-to-noise ratio V/N somewhat arbitrarily 
taken to be 3. Hence, if the ratio of photosensor output voltage (eq. (16)) to noise (eq. (21)) 
is set equal to 3 and solved for p(X), the minimum detectable albedo difference becomes 


<Ap) = 


3N 


— ^^D^kgRjG ^ S(X) ^c(^) ^(^) 


(46) 


MARS RADIANCE PROPERTIES 
Spectral Radiance 

The distance of Mars from the Sun varies from a maximum of 1.64 AU at the begin- 
ning of the mission to a minimum of 1.46 AU at the end of the mission (ref. 17). The 
solar irradiance at a distance of 1.6 AU is plotted in figure 15(a). 
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The atmosphere is considered to be a nonabsorbing "Rayleigh" atmosphere with nor- 
mal optical thickness T(X) plotted in figure 15(b) (ref. 17). The atmospheric transmit- 
tance is related to the optical thickness by 


r 

- 

- 


< -T('X) 

ij 

+ hg\ _ 2 ^ - cos t 

1 

ha 


X j 

/ 



- 

J 


(47) 


where rjjj is the planet radius and hg^ the equivalent atmospheric height of uniform 
optical thickness; the Mars radius rj^^ is 3400 km, and the equivalent atmospheric height 
hg is 25 km (ref. 18), Values of are plotted in figure 15(c) for several solar 

incidence angles 

Figure 15(d) presents plots of normal albedos p(X) that are representative for 
various areas on Mars (ref. 17). Perhaps still higher albedos should be anticipated in the 
blue-green wavelength range, since, according to a recent report (ref. 19), Martian color 
data obtained by the U.S.S.R. Mars 4 and 5 spacecraft indicate blue mountains and bluish- 
green craters. 


Illumination Scattering Fimction 

An illumination scatterii^ function formulated by Meador and Weaver (ref. 20) is 
used here to describe the diffuse reflectance of surfaces as a function of light incidence, 
emittance, and phase angles. The Meador-Weaver formulation is also a function of three 
parameters which are intended to represent physical surface properties (particle size, 
single-particle albedo, and compactness). However, no such significance is attached here 
to these parameters. They are used only as a convenient means of curve-fitting several 
illumination scattering characteristics with a single analytical function. Any dependence 
that the illumination scattering function can be expected to have on wavelength because of 
diffraction effects is not considered. 

The illumination scattering function is given by the expression 


0(t,e,g) = 


COS L 


+ aQ cos g^ f^t,e,g;a 2 ^ + aj(cos t + cos e) 


^1 + aQ + aj^(cos L + cos e) 
where f^t,e,g;a 2 ^ is a shadowing correction factor given by 


(48) 
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(t,e,g;a2) = + v ^ 

1 J 

exp< 

0 

6it 

3irx + 2(2 + x2j| 

1 ]l 

1 - x^ 1 + 6x sin“l X , 



L 


J 


dx 


Results from this function should be used only for values f^t,e,g;a2j = 1; all other values 
should be replaced by f^t,e,g;a2] = 1. The parameters fx and u are given by 


M = 


4a2(l + cos g) 
3 sin g 


7ra2(cos i + cos e)r » -|l/2 

i’ = sin^ g + 2(1 + cos g) cos t cos e 

sin g cos L cos e 

The latest and apparently best data on the illumination scattering characteristics of 
Mars are given by Thorpe (ref. 21) from Mariner 9 observations. Comparisons of these 
data with the Meador-Weaver function led to Uq = 0.55, a^ = 0.60, and a2 = 0.18 
(ref. 22); with these values the Meador-Weaver function predicts scattering very close to 
that given by Thorpe. Consequently, this set of parameters will serve here as nominal 
values. Three other sets of parameters are considered, as summarized in table V, pro- 
viding scattering characteristics that range from strong backscatter to near -Lambertian 
reflectance. 

Figures 16 to 18 present plots of the illumination scattering function for a wide 
range of lighting geometries. Figure 16 presents variations of the Mars (Thorpe) illumi- 
nation scattering characteristics with incidence angle for several emittance angles: 
€q = 40°, 50°, 60°, and 70°. The corresponding camera elevation aisles are -50°, -40°, 
-30°, and -20°, respectively. As illustrated in figure 2, these elevation angles encom- 
pass the surface sampler area. Figures 17 and 18 present variations of illumination 
scattering of the four different surfaces with incidence angle Iq and azimuth angle 9 
(see fig. 13), respectively. The emission angle €q is 60°, correspondir^ to a camera 
elevation angle of -30°, which points near the center of the sampler area. 

IMAGE QUALITY PREDICTIONS 
Constraints 

Predictions are presented to illustrate variations of image quality with lighting 
geometry for a wide range of surface radiance properties. The lander orientation with 
respect to the diurnal Sun path over the landing site imposes important constraints on 


27 



lighting and viewing geometries that are not accounted for by the predictions. Thus, 
otherwise favorable geometries may have to be avoided because of obscuration either by 
the lander structure or its shadow. 

The following constraints minimize the number of predictions necessary to cover a 
wide range of conditions: 

(1) The camera viewing geometry is constrained to an emittance angle Bq of 60°, 
or camera elevation angle of -30°, which points near the center of the sampler area (see 
fig. 2). 

(2) The selection of camera gains is constrained to those dynamic ranges that 
encompass the assumed surface radiance ranges. Table VI lists the maximum 
photosensor -array output voltages V for the four albedos shown in figure 15(d). 

(3) Broadband image quality predictions are presented as a function of the ratio 
V/N. Table VII presents these ratios for those combinations of surface albedos and 
camera gains that encompass the corresponding maximum signal levels. The value 
V/N = 50 is regarded to represent low signal-to-noise ratios; V/N = 100, average 
signal-to -noise ratios; and V/N = 200, high signal-to-noise ratios. 

(4) Cone diameters are expressed as a function of number of pixels c so that pre- 
dictions can be presented independently of object distance and instantaneous field of view. 
This simplification is only approximately valid. Defocus blur, for example, is a function 
of object distance. Also, the lens frequency response affects the smaller spatial detail 
associated with the narrow instantaneous field of view (/3 = 0.04°) more than the larger 
spatial detail associated with the broad instantaneous field of view (^ = 0.12°). 

(5) Predictions of cone signal-to-noise ratio (S/N) are constrained to cone diame- 
ters of 6 pixels and cone slopes of 10°, 30°, and 60°. Predictions of minimum detectable 
cone slopes (a) are constrained to cone diameters of 6 pixels. Predictions of minimum 
detectable cone diameters (c) are constrained to cone slopes of 30°. 

Broadband Image Quality 

Table VIII summarizes the broadband image quality predictions that are presented. 
Appendix C compares some of these predictions with experimental results. 

Figures 19 and 20 present graphs of the variation of cone signal-to-noise ratio 
(S/N), minimum detectable cone diameter (c), and minimum detectable cone slope (a) 
for the Mars (Thorpe) illumination scattering characteristics. Figure 19 presents these 
variations as a function of light incidence angle Lq for three azimuth angles {6 = 0°, 90°, 
and 180°). Figure 20 presents contour plots of (S/N), (c), and (o') as a function of inci- 
dence and azimuth angles. 
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Results show that the three broadband image quality criteria lead to generally sim- 
ilar conclusions with regard to favorable and unfavorable lighting geometry. However, 
figure 19 shows that (S/N) tends to stress favorable lighting near low Sim elevations, 
whereas both (c> and (o!> tend to stress unfavorable lighting near small phase angles 
(i.e., when the Sun is behind the camera, or 6 = 0°). 

Figures 21, 22, and 23 present graphs of the variation of cone signal-to-noise ratio 
(S/N) with lighting geometry for four illumination scattering characteristics. These var- 
iations are presented in figure 21 as a function of incidence angle Lq for three azimuth 
angles (6 = 0°, 90°, and 180°) and in figure 22 as a function of azimuth angle 0 for three 
incidence angles (lq = 40°, 60°, and 80°). Contour plots of (S/N) as a function of inci- 
dence and azimuth angles are presented in figure 23. 

Results suggest the following conclusions about image quality and lighting geometry: 

(1) Image quality depends strongly on the selection of a proper gain (i.e., dynamic 
range) for a given surface albedo. The dynamic range could be too low to encompass the 
scene radiance or too high to encode the scene radiance with an adequate number of quan- 
tization levels. Initial trials must be made and possible errors in selecting suitable gains 
must be anticipated. 

(2) Image quality depends strongly on lighting geometry. Sun elevation angles raid- 

ing from 10° to 30° are most often favorable. Lower Sun elevation angles, raiding from 
5° to 20°, tend to be more favorable for detecting shallow slopes, whereas higher Sun 
elevation angles, ranging from 20° to 40°, tend to be more favorable for imagii^ steep 
slopes and large details, such as rocks. Image quality diminishes more rapidly toward 
elevation angles lower than the most favorable ones rather than toward higher angles. It 
follows that it is better to err with Sun elevation angles that may prove higher than most 
favorable rather than lower. Shadows are then also less likely to obscure surface fea- 
tures. Important exceptions exist if the Sun is either behind the camera (i.e., 6 ~ 0°) or 

in front of it (i.e., 0 ~ 180°). 

If the Sun is behind the camera, then image quality can drastically diminish at small 
phase angles. Also important, although not a result of the predictions, is that Sun eleva- 
tion angles below the zero phase angle cause the target to be shadowed by the camera (and 
lander) structure. 

If the Sun is in front of the camera, then favorable Sun elevation angles are dependent 
on the target, low Sun elevation angles being preferable for cones with low slopes and high 
Sun angles for cones with steep slopes. It is also important to note, although not revealed 
by the curves, that the part of the target that faces the camera tends to be in shadow. It 
follows therefore that side lighting (40° % 9 % 120°) is favorable. 
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(3) Favorable lighting geometries are independent of surface albedo and do not vary 
appreciably for illumination scattering characteristics ranging from near-Lambertian 
reflectances to strong backscatter. However, surfaces with near-Lambertian reflec- 
tances provide a wider range of adequate lighting geometries than surfaces with strong 
backscatter. 


Narrowband Image Quality 

The variation of minimum detectable albedo difference (Ap) with lighting geometry 
is essentially dependent only on the surface illumination scattering function. The addi- 
tional dependence of this variation on atmospheric transmittance is small for Sun eleva- 
tion angles above 10°. Consequently, plots of the illumination scattering function pre- 
sented in figures 16 to 18 already reveal much about favorable lighting geometries. More 
complete information is presented in figure 24 by contour plots of (Ap) as a function of 
Sun azimuth and elevation angles for four illumination scattering characteristics. Numer- 
ical results in figure 24 are given for the green narrowband channel and a gaih number of 
4. This gain provides a dynamic range that encompasses all but the highest albedo shown 
in figure 15. Since the sensitivity of all narrowband channels is essentially limited by 
quantization noise for a gain number of 4, the numerical results approximately represent 
all narrowband channels. Furthermore, these values multiplied by 2 would approximately 
represent results for a gain number of 5. 

The curves show that high Sun elevation angles are always favorable. However, the 
preference shown for low phase angles as a result of backscatter is not generally valid. 

The construction of spectral reflectance curves from multispectral data may be adversely 
affected by nearly specularlike reflectances, and color pictures with very low surface 
contrast of both spatial and spectral detail may have a bleached appearance. 

CONCLUDING REMARKS 

Formulations were presented that permit predictions of image quality as a function 
of camera performance, surface radiance properties, and lighting and viewing geometry. 
The general objective of broadband (black-and-white) imaging to resolve small spatial 
details and slopes was formulated as the detectability of a right -circular cone with sur- 
face properties of the terrain. The general objective of narrowband (color and near- 
infrared multispectral) imaging to resolve spectral reflectance features was formulated 
as the minimum detectable albedo variation. The camera response was also formulated 
as a function of reference test chart features to aid in diagnosing its performance. A 
degraded performance is both diagnosed and simulated by changing camera response 
characteristics until image quality predictions agree with pictures of the reference test 
chart. 
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To satisfy primary objectives of the imaging investigation, it is important that the 
lander be so oriented that it avoids shadowing the sampler area at favorable Sun eleva- 
tions. A general goal is to obtain lighting of this area from two sides (i.e., morning and 
afternoon) so that all features can be viewed at one time or another out of shadow. 

Results of the image quality predictions suggest the following additional guidelines for 
imaging: 

1. Image quality depends strongly on the proper camera gain command for a given 
surface albedo. Improper gain commands could result in a dynamic range that is either 
too narrow to encompass the surface radiance or too wide to provide adequate quantiza- 
tion levels. 

2. Image quality depends strongly on lighting geometry. Side lighting with Sun ele- 
vation angles ranging from 10° to 30° is most often favorable for broadband imaging. 
Lower Sun elevation angles, ranging from 5° to 20°, tend to be more favorable for detect- 
ing shallow slopes, whereas higher Sun elevation angles, ranging from 20° to 40°, tend to 
be more favorable for imaging steep slopes and large details, such as rocks. High Sun 
elevation angles are always favorable for narrowband multispectral imaging. 

3. Favorable lighting geometries are independent of surface albedo and do not vary 
appreciably for a wide range of illumination scattering functions. However, surfaces with 
near-Lambertian reflectances tend to provide a wider range of adequate lighting geom- 
etries than surfaces with strong backscatter. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
March 2, 1976 
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APPENDIX A 


RADIOMETRIC RESPONSE 

The purpose of this appendix is to present a complete set of data required to predict 
camera radiometric response and sensitivity, and to describe absolute radiometric cam- 
era calibrations. Data are presented here for only one flight camera (serial no. 08); how- 
ever, the data are also approximately valid for the other flight cameras. A complete set 
of data for the four flight cameras is given in reference 23 together with their absolute 
radiometric calibration data. 

Table IX presents optical throughput and photosensor -array responsivities; table X 
gives photosensor -array electrical characteristics. The gain constant of the video pro- 
cessing electronics kg is 444.321 N(j/V, the offset interval k(.o is 0.1441 V, and the 
fixed offset ko is 0.204 V. 

Absolute radiometric calibrations were made by imaging a reference test chart 
which, in turn, was illuminated by a lamp that had been calibrated by Eppley Laboratory, 
Inc., the National Bureau of Standards (NBS). A special calibration fixture was used to 
insure that the lamp-to-chart distance and the lighting and viewing angles remained con- 
stant for all measurements. Four major error sources and the errors are as follows: 


Lamp irradiance ±3% 

Reference test chart reflectances .... ±6% 

Fixture ±3% 

Camera gains and offsets ±3% 

Root -sum-square error ±8% 


The lamp irradiance was calibrated to an accuracy of ±3 percent in the spectral range 0.4 
to 1.1 jiim (ref. 24). 

The reference test chart reflectances were measured relative to magnesium carbo- 
nate (MgC 03 ), for which the absolute reflectance in this spectral range is known with an 
accuracy of ±3 percent (ref. 25). Errors introduced by not carefully accounting for some 
of the variations of this reflectance with wavelength and with lighting and viewing geome- 
try, both in measuring the reference test chart reflectances and in using these reflectance 
data, diminish the accuracy to ±6 percent. To measure and use reference test chart re- 
flectances more accurately would significantly increase calibration complexity. The 
reflectances of the 11 gray patches are given in figure 12. These reflectances rep- 
resent the average values of measurements of three reference test charts; the measured 
reflectances varied less than ±3 percent from these average values (ref. 26). These vari- 
ations between charts are within measurement errors and are neglected. 
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APPENDIX A 


The calibration fixture was calibrated to account for peculiarities of lighting geome- 
try and for internal reflections which occur despite careful baffling. The remaining un- 
certainty was estimated to be ±3 percent. 

The gain and offset constants kg, k^o, and ko vary independently with gains, off- 
sets, and temperature (-410 to +lQO C) by less than ±2 percent around their average value 
for each camera; the total error is about ±3 percent. One exception exists: The gain 
constant kg for a gain number of 0 differs much more from this value; however, this 
gain was not used to obtain calibration data. 

Predictions made by using equations (13), (16), and (17) and calibration measure- 
ments agree if the calibration constants kg listed in table X are used in equation (13). 
These calibration constants reveal that the measured photosensor -array output voltages 
are consistently higher than their predicted values. These discrepancies apparently arise 
from differences between the methods used for calibrating the photosensor array and in 
calibrating the cameras, and from light reflections internal to the array (ref. 23). 
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APPENDIX B 


SPATIAL RESOLUTION 

The purpose of this appendix is to compare predictions and measurements of the 
image contrast of three tribars located on the reference test charts. The capability to 
predict image contrast of known spatial detail aids in diagnosing camera performance. 

Figure 25 shows three Viking lander camera images of the reference test chart. 

The images were obtained with the two nearest -focused high -resolution diodes (BBl and 
BB2) and the survey diode. The camera was located 1.0 m from the chart and viewed it 
normally. Figure 26 compares the normalized contrast of these tribar images with the 
normalized camera square -wave response AV/V(t) predicted by equation (23). Predic- 
tions and measurements are in good agreement, the larger tribar being distinctly repro- 
duced in all three images. Contrast of the medium tribar is slightly reduced in the BBl 
image, significantly reduced in the BB2 image, and reduced to near threshold in the survey 
image. The smallest tribar is resolved at a very low contrast in the BBl image, and gives 
rise to a so-called false resolution in the other two modes (i.e., in this case two rather 
than three bars). 


^4 
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APPENDIX C 


BROADBAND IMAGES 

The purpose of this appendix is to compare predictions and nr^easurements of the 
cone signal -to -noise ratio. Results are intended to verify predicted variations of image 
quality with lighting geometry, and to illustrate some of the limitations that are inherent 
in the use of a single criterion to describe the quality of an image. 

Figure 27 presents comparisons of predicted and measured cone signal -to -noise 
ratios (S/N), and images from which the measurements were obtained; figure 27(a) is for 
a Sun azimuth angle 9 = 40° (which was the closest to 0=0° that could be obtained for 
measurements), figure 27(b) is for 9 = 90°, and figure 27(c) is for 9 = 180°. Measure- 
ments of (S/N) were obtained at 10° intervals of incidence angle Lq ranging from 30° to 
80°. A cone with slope a = 30° and radius rcn = 3 cm was located at a distance of 
2.7 m, corresponding to the in-focus distance of the high -resolution diode BB2. The sur- 
face material was Colorado basalt (latite), with a mean particle diameter of 220 p.m. The 
target was illuminated by sunlight with the aid of a heliostat. (The curves shown in fig. 27 
depart slightly from corresponding curves shown in fig. 19 at high incidence angles be- 
cause atmospheric transmittance is independent of incidence angle when the heliostat is 
used.) 

Evaluations of experimental results were based on computer printouts of the image 
data that gave the quantization level of each pixel. The evaluator first outlined the part of 
the cone that is brighter than the background and the other part that is darker (including 
shadow, if present). Then he recorded the number of pixels n± in each section and their 
quantization level N^^i. The effective diameter of each cone section was taken to be 



and the average signal excursion was taken to be 


n± 

I |Nd,l-Nd.o| 

i=l 

where Njj^q is the baci^round quantization level. The corresponding photosensor-array 
output voltages V± were computed with equation (17); the signal-to-noise ratios (S/N), 
with equation (45a). 
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Results show that measurements are in general agreement with Image quality pre- 
dictions. However, the images also illustrate an obvious limitation of the predictions in 
failing to account for the orientation of surface target shadows: Although high surface 
contrasts can be obtained when the camera faces the Sxm, shadows obscure the side of the 
cone (i.e., rock) that could otherwise be observed. 
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TABLE I.- DESIGN AND OPERATING CHARACTERISTICS OF CAMERA 


Height, cm 

Diameter: 

Upper elevation assembly, cm 

Lower azimuth assembly, cm 

Weight, kg 

Power: 

Internal calibration, W 

Low-resolution imaging, W 

High-resolution imaging, W 

Standby, W 

Thermal control (av for cold day), W . . 
Slew to new command (less than 30 ms), W 

Environment: 

Atmosphere 

Atmospheric pressure: 

Earth, mbar 

Cruise 

Mars, mbar 

Ambient temperature: 

Sterilization, 

Mars, °C 

Hazards 


55.6 

14.4 

25.6 

7.26 

22 

34 

27 

22 

3 

111 


Earth or CO 2 

1000 

Vacuum 

2.8 to 20 

113 

-126 to +52 

Dust storms up to 70 m/s 



TABLE II.- PERFORMANCE CHARACTERISTICS 


Characteristic 

Survey 

Color and IR 

High resolution 

Instantaneous field of view, deg 

0.12 

0.12 

0.04 

Picture element registration error, deg 

±0.036 

±0.013 

±0.006 

Absolute angle error: 




Elevation, deg 

±0.3 

±0.2 

±0.2 

Azimuth, deg 

±0.15 

±0.1 

±0.1 

Frame width: 




Elevation, deg 

61.44 

61.44 

20.48 1 

Azimuth (min; max), deg 

2.5; 342.5 

2.5; 342.5 

2.5; 342.5 i 

Field of view: 



! 

1 

Elevation, deg 

100; from 40° above to 60° below horizon in 10° steps 

Azimuth, deg 

342.5; in multiples of 2.5° steps 1 

Geometric depth of field, m 

1.7 to 

1.7 to 

1.7 to °° 

In -focus distance, m 

3.7 

3.7 

1.9, 2.7, 4.5, and 13.3 

; Picture elements per line 

512 

512 

512 

j 

i Bits per picture element 

6 

6 

6 

' Bits per degree azimuth 

2.84 X 104 

8.53 X 104 

8.53 X 104 

i Time per degree azimuth: 




Rapid scan, s 

1.84 

5.52 

5.52 

Slow scan, min 

2.0 

6.0 

6.0 



TABLE in,- TYPICAL PICTURE FORMATS: THEIR NUMBER 


OF BITS AND TIMES NECESSARY TO OBTAIN AND 
TRANSMIT DATA AT 16 000 bits/s 

Average amoimt of imagery data transmitted per day: 
Nominal mission: 10 x 10® bits recorded, 6 x 10® bits 
real time 

Preprogramed mission: 10 x 10® bits 


Imaging mode 

Survey: 

60° X 60° 

60° X 120° 

60° X 240° 

60° X 360° 

Multispectral (color or IR) : 
60° X 60° 

60° X 120° 

High resolution: 

20° X 20° 

20° X 60° 

200 X 1200 


Time, min 

Data, bits 

1.84 

1.7 X 10® 

3.68 

3.4 

7.36 

6.8 

11.04 

10.2 

5.52 

5.1 X 10® 

11.04 

10.2 

1.84 

1.7 X 10® 

5.52 

5.1 

11.04 

10.2 



TABLE IV.- TYPICAL NOISE -EQUIVALENT RADIANCES 
FOR ALL GAINS AND SCAN RATES 


Scan 

rate 

Gain 

number 

Noise-equivalent radiance, w/m^-sr 

for channel - 


High 

resolution 

Survey 

Blue 

Green 

Red 

mi 

m2 

IR3 

Slow 

0 

0.013 

0.008 

0.014 

0.013 

0.010 

0.011 

0.013 

0.011 


1 

.019 

.016 

.020 

.019 

.018 

.018 

.019 

.017 


2 

.035 

.032 

.036 

.032 

.035 

.033 

.033 

.030 


3 

.067 

.064 

.068 

.062 

.070 

.065 

.063 

.059 


4 

.133 

.128 

.135 

.122 

.141 

.130 

.124 

.118 


5 

.265 

.256 

.270 

.243 

.281 

.259 

.247 

.235 

Rapid 

0 

0.069 

0.019 

0.080 

0.076 

0.040 

0.055 

0.073 

0.058 


1 

.070 

.024 

.081 

.077 

.043 

.056 

.074 

.059 


2 

.076 

.036 

.086 

.081 

.053 

.063 

.079 

.064 


3 

.095 

.066 

.104 

.097 

.080 

.084 

.095 

.082 


4 

.149 

.129 

.156 

.143 

.146 

.140 

.143 

.131 


5 

.274 

.257 

.281 

.254 

.284 

.264 

.257 

.242 




TABLE V.- ILLUMINATION SCATTERING PARAMETERS OF FOUR SURFACES 


Surface 

Physical parameters 

ao 

H 

^2 

Mars (Thorpe) 

0.55 

0.60 

0,18 

Colorado basalt (mean particle diameter of 220 ju.m) 

-.10 

.26 

.15 

Strong backscatter 

.8 

.1 

.1 

Near -Lambertian reflectance 

-.8 

3.0 

i 

.5 


TABLE VI,- MAXIMUM SIGNAL LEVELS FOR FOUR SURFACE ALBEDOS 


Maximum photosensor -array output voltage, V, for channel — 


Surface 

albedo 

Broadband 


Color 


Near -infrared 


Blue 

Green 

Red 

IRl 

IR2 

IR3 

Pi 

2.40 

1,87 

2.45 

2.31 

2.42 

2.57 

2.72 

P2 

1.74 

1.43 

1.81 

1.68 

1.82 

1.93 

2.03 

P3 

1.24 

1.24 

1.39 

1.15 

1.27 

1.35 

1.40 

P4 

.89 

1.11 

1.14 

.86 

.77 

.81 

.89 


TABLE Vn.- V/N FOR SEVERAL SURFACE ALBEDOS 
AND BROADBAND IMAGING MODES 


Gain number 

Imaging mode 

V/N for albedo - 

Pi 

P2 

P3 

P4 

3 

High resolution, rapid scan 




121 


High resolution, slow scan 




172 


Survey, rapid scan 




173 

4 

High resolution, rapid scan 


151 

108 

77 


High resolution, slow scan 


170 

121 

87 


Survey, rapid scan 


175 

124 

89 

5 

High resolution, rapid scan 

111 

82 

59 

42 


High resolution, slow scan 

115 

85 

60 

43 


Survey, rapid scan 

119 

88 

62 

45 


44 


TABLE Vin.- SUMMARY OF BROADBAND IMAGE QUALITY PREDICTIONS 


Figure 

Image -quality 
criteria 

V/N 

Lighting geometry^ 

Cone 

Scattering 

characteristics*3 

Incidence, 
to , deg 

Azimuth, 
0, deg 

Slope, 
a, deg 

Diameter, 
c, pixels 

19 

(S/N) 


Variable 

0, 90, 180 

10, 30, 60 

6 

Mars (Thorpe) 


<c) 

50, 100, 200 



30 

— 



<«> 

50, 100, 200 




6 


20 

<S/N> 

100 

Variable 

Variable 

30 

6 

Mars (Thorpe) 


(c> 




30 

— 



(a) 





6 


21 

(S/N> 


Variable 

0, 90, 180 

10, 30, 60 

6 

All four 

22 

(S/N) 


40, 60, 80 

Variable 

10, 30, 60 

6 

All four 

23 

<S/N) 


Variable 

Variable 

30 

6 

All four 


^See figure 13. The emittance angle eo is 60°. 
*^See figures 16 and 17. 
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TABLE IX ,- TRANSMITTANCE OF OPTICS AND RESPONSIVITIES OF PHOTOSENSOR ARRAY 


Transmittance Responsivity , A / W , of channel 


I 

i 1 

^ cw > 

y 


BBl 

BB2 

BB3 

BB4 

Survey 

Blue 

Green 

Red 

IRl 

IR2 

IR3 

0.400 

0,926 

0.692 

0.931 

0.096 

0.095 

0.096 i 

0.104 

0.097 

0.044 

0.001 

0.005 

0.002 

0.002 1 

0.005 

.425 

.930 

.774 

,945 

.127 

.125 

.129 

.132 

.126 

.082 

.001 

.004 

.001 

.001 

.003 

.450 

,932 

.783 

.959 

.150 

.149 

,159 

.163 

.151 

.120 

.001 

.002 

.001 

.001 

.005 

.475 

.934 

.820 

.964 

.174 

.172 

.187 

.188 

.175 

.142 

.003 

.001 

.001 

.002 

.009 

.500 

.940 

' .857 

.963 

.197 

.196 

.212 

.215 

.198 

.120 

.044 

■ .001 

.000 

.008 

1 .014 

.525 

.940 

.846 

.966 

.204 

.201 

.221 

.224 

.210 

.004 

.108 

.001 

.000 

.017 

.008 

.550 

.945 

.844 

.954 

.214 

.215 

.235 

.240 

.224 

.001 

.186 

.001 

.000 

.007 

.005 

.575 

.945 

.837 

.953 

.232 

.236 

.259 

.266 

.248 

.001 

.082 

.003 

.001 

.001 

.005 

.600 

.949 

.817 

.955 

.247 

.253 

.279 

.286 

.267 

.001 

.017 

.053 

.001 

.001 

.006 

.625 

.948 

.819 

.950 

.264 

.270 

' .297 

. .306 

.285 

.002 

.003 

.252 

.001 

.001 

.008 

.650 

.950 

.800 

.943 

.276 

.282 

.310 

.317 

.296 

.014 

.001 

.251 

.001 

.001 

.010 

.675 

.949 

.795 

.942 

.282 

.290 

.321 

.327 

.306 

.037 

.001 

.265 

.001 

.001 

.015 

.700 

.949 

.790 

.947 

.291 

.298 

.330 

.335 

.316 

.005 

.001 

.293 

.000 

.001 

.022 

.725 

.948 

.767 

.934 

.297 

.307 

.338 

.345 

.326 

.001 

.001 

.187 

.001 

.002 

.029 

.750 

.948 

.744 

.945 

.307 

.315 

.345 

.352 

.332 

.001 

.000 

.040 

.001 

.005 

.042 

.775 

.947 

.749 

.926 

.315 

.323 

.354 

.364 

1 .343 

.000 

.002 

.024 

.006 

.007 

.007 

.800 

.948 

.731 

.920 

.320 

.328 

.359 

.367 

.351 

.002 

.002 

.010 

.030 

.002 

.002 

.825 

.947 

.785 

.916 

.321 

.330 

.364 

.375 

.358 

.000 

.000 

.007 

.128 

.002 

.002 

.850 

.947 

.782 

.914 

.319 

.327 

.358 

.370 

.356 

.000 

.000 

.003 

.290 

.002 

.002 

.875 

.944 

.794 

.904 

.306 

.314 

.343 

.353 

.344 

.000 

.000 

.007 

.226 

.005 

.003 

.900 

.943 

.815 

.895 

.286 

.293 

.319 

.329 

.328 

.001 

.000 

.001 

.272 

.025 

.008 

.925 

.945 

.835 

.892 

.265 

.272 

.297 

.305 

.309 

.000 

.001 

.000 

.101 

.220 

.028 

.950 

.945 

.849 

.890 

.240 

.244 

.267 

.275 

.285 

.001 

.001 

.000 

.026 

.270 

.117 

.975 

.947 

.863 

.882 

.209 

.213 

.235 

.239 

.253 

.003 

.001 

.000 

.011 

.138 

.212 

1.000 

.941 

.870 

.869 

.159 

.162 

.179 

.181 

.195 

.002 

.003 

.000 

.006 

.027 

.169 

1.025 

.945 

.876 

.873 

.112 

.114 

.128 

.127 

.140 

.002 

.002 

.001 

.004 

.010 

.121 

1.050 

.942 

.880 

.862 

.067 

.067 

.078 

.075 

.083 

.004 

.002 

. .001 

.003 

.005 

.072 

1.075 

.942 

, .880 

.861 

.035 

.034 

.040 

.038 

.043 

.018 

.006 

.000 
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TABLE X.- ELECTRICAL CHARACTERISTICS OF PHOTOSENSOR ARRAY 


Channel 

kc 

Rf, MS2 

Ri, 

Channel 

gain 

rp, iim 

Diode noise, 
fA/\/Hz 

Preamplifier noise, 
fA/\/Hz 

Total noise, 
fA/\/Hz 

BBl 

1.18 

728.5 

5 

178 

19.81 

19.4 

5.98 

4.97 

7.78 

BB2 

1.17 

699.7 

5 

680 

17.61 

20.8 

6.25 

5.09 

8.07 

BBS 

1.18 

724.9 

5 

958 

16.78 

19.7 

6.06 

4.80 

7.74 

BB4 

1.19 

740.8 

6 

480 

15.48 

19.4 

5.98 

5.25 

7.97 

Survey 

1.18 

756.8 

54 

540 

1.88 

58.8 

18.09 

4.96 

18.76 

Blue 

1.85 

785.0 

4 

892 

22.77 

59.4 

4.97 

4.58 

6.77 

Green 

1.89 

709.5 

4 

082 

24.50 

58.8 

4.64 

5.07 

6.89 

Red 

1.07 

764.1 

14 

240 

7.02 

59.1 

8.77 

4,81 

10.01 

IRl 

1.09 

719.5 

7 

064 

14.16 

58.2 

6.14 

4,98 

7.92 

IR2 

1.26 

707.1 

4 

561 

21.98 

58.8 

5.18 

5.17 

7.88 

IRS 

1.82 1 

752.2 

6 

087 

16.56 

58.8 

5.75 

4.88 

7.52 


<1 


4 ^ 

00 
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Figure 1.- Artist’s view of Viking lander. 
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Figure 3.- Simplified cutaway view of Viking lander camera. 




Figure 4.- Elevation scan profile. Angular and time values are given for the high- 
resolution (X = 0.04°) and rapid-scan (4.7 scans/s) mode; angular values are to be 
multiplied by 3 for low -resolution (X = 0.12°) imaging, and time values are to be mul 
tiplied by 64 for the slow scan rate (0.073 scans/s). 
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Figure 5.- Spectral responsivities of photosensor array. 




•PHOTOSENSOR ARRAY VIDEO PROCESSING ELECTRONICS- 


Figure 6.- Simplified circuit diagram of photosensor array and video electronics 







Angular frequency Jp/deg 


Figure 9.- Spatial frequency response, u = 45 represents maximum defocus blur for 
low-resolution imaging (/3 = 0.12°) over the range from 1.7 m to u = 15, for 
high-resolution imagii^ 0 = 0.04°) over the same range if the best-focused diode 
is used at each distance. 



(c) Lens. 



(d) Camera. 
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Figure 10.- Optical throughput. 
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Figure 12.- Reference test chart. (Numbers 
in patches are reflectances.) 
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(a) Geometry for incident (b) Geometry for emitted (c) Top view of cone geometry 

radiation. radiation. 

Figure 14.- Cone geometry. 
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Figure 19,- Variation of image quality ((S/N>, (c), and (a>) with light incidence angle t 
for Thorpe’s illumination scattering function of Mars. Azimuth angles 6 are 0°, 
90°, and 180°; emittance angle eo is 60° (i.e., camera elevation angle of -30°), 




(a) Cone signal-to-noise ratio (S/N>. a = 30°; c = 6 pixels. 



(b) Minimum detectable cone diameter (c>. a = 30°. 



(c) Minimum detectable cone slope (a), c = 6. 

Figure 20.- Contour plots of image quality as a function of light incidence angle Iq (or 
Sun elevation angle SEA) and azimuth angle 6 for Thorpe’s illumination scattering 
function of Mars. Emittance angle 6o is 60° (i.e., camera elevation angle of -30°); 
imaging condition is adequate, V/N = 100. 
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Figure 21.- Variation of cone signal-to-noise ratio (S/N) with light incidence angle Lq, Plots are given 
for three azimuth angles 6 and four illumination scattering characteristics. 






(a) Mars (Thorpe). 


(c) Colorado basalt. 



(b) Strong backscatter. (d) Near Lambertian. 


Figure 23.- Contour plots of cone signal-to-noise ratio (S/N) as a function of light incidence angle Lq 
(or Sun elevation angle SEA) and azimuth angle 6 for four illumination scatterii^ characteristics. 
Emittance angle 6 q is 60® (i-Gv camera elevation angle of -30®); cone slope o! is 30®; diameter 
c is 6 pixels; imaging condition is adequate, V/N = 100. 
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(d) Near Lambertian. 


Figure 24.- Contour plots of minimum detectable albedo difference (Ap) as a function of light incidence 
ai^le Lq (or Sun elevation angle SEA) and azimuth angle 9, for four illumination scattering char- 
acteristics. Emittance angle Cq is 60° (i.e., camera elevation angle of -30°). Values are for a 
gain number of 4; twice these values represent a gain number of 5. 
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Figure 25.- Images of the reference test chart obtained with the two nearest-focused 

high-resolution diodes and the survey diode. 
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Figure 26.- Predicted square-wave frequency response of the camera for a 
target located 1.0 m away, and normalized tribar contrast measurements 
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(a) e = 40°. 


Figure 27.- Cone images and their predicted and 
measured signal -to -noise ratios. 
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(c) e = 180°. 

Figure 27.- Concluded. 
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